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' Abstract 



We introduce a new approach to highly correlated systems which general- 
izes the Fermi Hypernetted Chain and Correlated Basis Function techniques. 
While the latter approaches can only be applied to systems for which a non- 
relativistic wave function can be defined, the new approach is based on the 
variation of a trial hamiltonian within a path integral framework and thus 
can also be applied to relativistic and field theoretical problems. We derive a 
diagrammatic scheme for the new approach and show how a particular choice 
of the trial hamiltonian corresponds exactly to the use of a Jastrow correlated 
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ansatz for the wave function in the Fermi Hypernetted Chain approach. We 
show how our new approach can be used to find upper bounds to ground state 
energies in systems which the FHNC cannot handle, including those described 
by an energy-dependent effective hamiltonian. We demonstrate our approach 
by applying it to a quantum field theoretical system of interacting pions and 
nucleons. 

21.60.-n;11.10.Ef;21.65.+f 
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I. INTRODUCTION 



Variational methods provide extremely powerful tools for analyzing complicated many- 
body hamiltonians. Fermi Hypernetted Chain (FHNC) techniques, in particular, have been 
successfully applied to many different problems, ranging from nuclear matter to strongly 
interacting electrons to liquid helium, where strong correlations between particles dominate 
the system and prohibit the usual perturbative expansions0-i. Correlated Basis Function 
(CBF) theory, which employs FHNC and cluster summation methods to compute pertur- 
bative corrections to quantities determined variationally, has also been successfully applied 
to a wide variety of systems@-i. These methods, however, have up to now been limited 
to hamiltonian systems for which a nonrelativistic wave function can be defined. In the 
present work we build upon a previous letteii, developing a framework which goes beyond 
this limitation. Our particular goal is to find a variational method for studying nuclear 
systems which include dynamical meson exchange. The methods we shall develop, however, 
are much more general. 

The basic idea of our approach is to adopt a Feynman path integral formulation in 
euclidean space so as to choose a trial hamiltonian which may be determined variationally. 
This is less restrictive than the usual FHNC approach, where one chooses a variational wave 
function. The use of a variational principle in Feynman path integrals is nothing new, as we 
describe in Sec. 2. The novelty is connecting this principle to the extensive existing FHNC 
calculations. This is a major goal of the present work. We shall show that the FHNC 
diagrammatic expansion can be recovered as a particular case of a more general expansion 
based upon Feynman diagrams. Once this contact is made, one may develop extentions 
of the highly successful FHNC methods for summing up diagrams to be applied in the 
more general path integral context, thus opening an aspect of the path integral variational 
method little explored up to now. Such techniques should have a wide range of applications, 
including possible improvements upon many existing FHNC calculations, as well as the 
analysis of systems that the usual FHNC cannot handle, such as those described by effective 
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(time-dependent) hamiltonians. 

The path integral formulation has another advantage: it yields a perturbative expansion 
which is more compact, as well as more general, than that of time-independent perturbation 
theory. This can be a great technical advantage — for example, one Feynman diagram with 
n interaction lines corresponds in principle to n! different Goldstone diagrams. The path 
integral formulation thus provides a convenient setting for the study of strongly interacting 
systems: the variational principle can be used to find a trial hamiltonian which reproduces as 
closely as possible the physics of the true hamiltonian, and the corrections can be calculated 
within a natural perturbative framework. Since it is usually a matter of some art to choose 
a trial hamiltonian which successfully balances the competing demands of accuracy and 
simplicity, a framework which allows a ready interplay between the determination of a 
variational minimum and perturbative corrections is clearly very useful. 

The essential simplification which we shall adopt in choosing a trial hamiltonian is mo- 
tivated by the FHNC expansion. Whereas a normal perturbative approach builds collective 
states from independent particles, the FHNC instead treats interparticle correlations as fun- 
damental. Our idea is thus to choose a trial hamiltonian hg containing a nontrivial two-body 
piece which may be varied to reproduce as best as possible the interactions in the true hamil- 
tonian h, but we shall require the one-body part of ho to have a highly simplified form with 
respect to that of h. In contrast to this, the usual perturbative expansion takes ho to be just 
the one-body kinetic energy operator, with no two-body term. Our choice, instead, gives 
the interactions a more fundamental role; and the simplification of the one-body part shall 
allow us, as in the FHNC expansion, to sum entire classes of diagrams which give important 
contributions in strongly correlated systems. 

As a matter of fact, we shall choose a one-body piece which depends on a single parameter 
a, and we shall demonstrate that the choice a — (which corresponds to "turning off" 
the one-body piece completely) exactly reproduces the FHNC expansion. So one immediate 
result of our formulation is that we can now calculate upper bounds on effective hamiltonians 
by using the full machinery of the usual FHNC approach. The choice a = is essentially a 
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static approximation, and to go beyond this approximation one must either take a > or 
include higher order corrections in the perturbation h — ho. If we take a > 0, we must extend 
the FHNC techniques for summing diagrams, because now the time integrations are no longer 
trivial. While it is surely possible to use finite temperature techniques to generalize the 
FHNC diagrammatic summation so that a may be treated as a true variational parameter, 
in the present work we shall be content with looking at the particular case a oo. In this 
limit, the time integrations simplify considerably, and we shall show a particular example of 
how to sum up diagrams in analogy with the FHNC expansion. 

For the particular case a = 0, the trial hamiltonian ho is directly related to a trial wave 
function of the CBF type, with the correlations in the CBF wave function corresponding 
directly to the interactions in Hq. Of course, as with the CBF correlations, there is a lot 
of freedom in choosing the two-body part of /iq- Obviously, the closer one can make this 
to that of h, the better one's variational estimate of the energy will be. The inclusion 
of state dependence (that is, spin and isospin dependence) in Hq clearly improves varia- 
tional estimates, although in the usual FHNC there are nevertheless many difficulties in 
implementing this generalization. A popular approximation is the so-called Single Opera- 
tor Chain approximationi, which consists of summing chains of single operator correlations 
without hyperconnecting them. However, it seems necessary to go beyond the SOcil, and it 
is not yet well understood how to do this. A new approach using spin coherent-state wave 
functions currently being developed0 is one possibilty. While the present approach may very 
well provide an alternative means of improvement in this direction, either by allowing for a 
more complete summation of state-dependent correlations through techniques developed for 
nonzero a or by allowing for a more straightforward calculation of state-dependent correc- 
tions through perturbation theory, we shall limit ourselves in the following to an Hq which 
contains only spin-independent two-body interactions. This corresponds to using the the 
Jastrow (state-independent) correlated trial wave function ansatz. The study of how best 
to include state-dependent correlations will thus be an important immediate application of 
the approach we develop here. 
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The outline of the paper is as follows. We describe the variational principle in the 
path integral setting and discuss the choice of the trial hamiltonian in Sections 2 and 3. 
Sections 4 through 6 are concerned with obtaining the usual FHNC diagrammatic scheme 
as a particular limit of a more general Feynman diagrammatic scheme. In the process of 
demonstrating this, we show that simplifications which allow the use of FHNC techniques 
for summing diagrams do not occur in general. In Section 7, in a different limit from that 
which recovers the FHNC diagrammatic scheme, we show how diagrams may be summed in 
a different way, and hint at generalizations of the FHNC techniques for this case. In Section 
8 we demonstrate our procedure by applying it to a nuclear system containing dynamical 
pions. Finally, in Section 9 we discuss our results and point toward directions for future 
study. 

II. THE VARIATIONAL PRINCIPLE IN TERMS OF FEYNMAN PATH 

INTEGRALS 

The original idea of formulating the variational principle in the path integral context is 
due to Feynman himselfl'lli. Let us briefly summarize here the formalism. 

As is well known, the ground state energy of a system can be expressed in terms of the 
partition function according to 

E„= ii,„!^. (1) 

On the other hand, the partition function can be expressed in a path integral form as 



where we have assumed that the system can be described in terms of a field ip (either 
fermionic or bosonic) and its hermitian conjugate. Here h[ip'^ , ip] is the classical hamiltonian. 
More precisely, when some ambiguities arise in performing the Legendre transformation, we 
should start from the lagrangian of the system L[^'^{t)^ip{ty\ and make the transformation 
to the euclidean world (imaginary time) t ^ it = t. 
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Once a representation of Z is provided, we can define an average over tlie functional 
space where the Feynman-Kag integral is defined, namely, 

( / ).o = ^/ /^[v^^^]e-/-^>°t^^*7[^^^] , (3) 

with 



The weight function e ^-i^i'^^^'^^ ''^''^^ must be positive definite, which only implies that /iq 
must be real. Using this definition the partition function can be rewritten in the form 



\ I ho 

Now, it is a general property of a mean that 



e^) > e<^^ . (6) 



Thus, applying this last inequality to and inserting into ([^), we find 
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-log! /D[^t,^]e-/-^>o'^-) (7) 



< lim 

v-^/'' /ho y-^i^ I ho. 

The first term on the r.h.s. of (|^) cancels with the third. This outcome is rather obvious 
if /lo is time independent. In this case the first term is by definition the ground state 
energy of the hamiltonian (being time independent and real, /iq can in fact be regarded 
as a hamiltonian) and the third term is the expectation value of the same hamiltonian in 
the ground state, hence the cancellation. However, since we are interested in extending the 
variational scheme to effective hamiltonians and therefore want to maintain the possibility 
of a time dependence in /iq, we need a more refined proof, which we present in Appendix A. 
There we show that even for time-dependent potentials, we end up with 

Eo< lim - / h{r)dr) . (8) 
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Reexamining the derivation of we see that the functional integral describing Z re- 
quires, in order to ensure its existence, that h must be real and bounded from below. A 
possible dependence upon the time in h (or conversely, in energy-momentum space, upon 
the energy) cannot be ruled out on the ground of mathematical reasons. We know that 
the hamiltonian is energy-independent indeed, but it can be useful sometimes to derive an 
effective hamiltonian, whose energy-dependence is the price we pay for having eliminated 
some unwanted degrees of freedom. Even in this case the previous procedure applies safely 
with the only constraint that the effective hamiltonian (obviously bounded from below be- 
cause the original one is) must be real in the euclidean world (that is, after the replacement 



The same consideration applies to as well, which is nothing but a weight function, even 
if it will be referred to in the following as a "trial hamiltonian" : again, the only condition 
/lo must fulfill for the variational principle (|]) to be mathematically well-defined is that it 
must be real. Of course, for the (5 oo limit to be physically meaningful, we require that 

must have an unique ground state. This is related to the choice of boundary conditions 
in (|^), as we shall see in Sec. 4. For the present we simply note that the inequality is valid 
independent of the boundary conditions on the path integral, and indeed also for finite (5. 
Thus we have a variational principle 



for any Z of the form (^, and the boundary conditions and /? — oo limit must be chosen 
so as to make Fq correspond to the quantity in which we are interested, namely, the ground 
state energy. 



in Feynman diagrams. Because of the exponentiation property of connected diagrams, we 
may write immediately 



E iE). 





Starting from (^, we may also find perturbative corrections in hj = h ~ ho to the 
variational energy (P) by simply expanding the ratio 




(10) 
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where (...)*" denotes connected diagrams. The quantity W = In Z = In Zq + In Z' is thus the 
sum of connected Feynman diagrams containing no hj interaction hues, one hj line, two hj 
hues, and so on. One easily sees that to first order — hm W/(3 is the variational energy (Rl), 
and the higher orders bring in the corrections. There will of course be a delicate interplay 
between how far one wants to proceed in optimizing the trial hamiltonian and how far one 
wants to be able to go in calculating perturbative corrections, decisions which will depend 
greatly upon the particular system one is studying and upon the particular quantities one is 
interested in calculating. We see here that the Feynman path integral formulation provides 
a convenient framework for arriving at such decisions. 



III. THE TRIAL HAMILTONIAN 

Different choices of Hq lead to different approximation schemes. For instance, we could 
choose ho in the form of a single-particle kinetic energy operator. It is immediately seen that 
in such a case the variational method coincides with the first order of the usual perturbative 
expansion. We could however assume ho in the form of a single particle operator but retain 
a possible time dependence. This topic could give rise to interesting investigations, because 
very much in the same line Feynman investigated the problem of the polaron, obtaining 
numerical results far better than the ones provided by the standard techniques (perturbation 
theory and canonical transformations)^. 

The problem can also be handled at a different level, by introducing a more complicated 
trial hamiltonian having not just a single particle but also a two particle term. The con- 
dition we require, of course, is that these hamiltonians can be handled, technically, much 
more easily than the full hamiltonian of the system. Since a many body system is intrin- 
sically complicated by the presence of the two-body interaction, no matter how simple this 
interaction is, the possibility we might explore instead is to simplify as far as possible the 
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single particle term. In this respect two options naturally arise, as we shall discuss in the 
following. 

Consider first a general lagrangian of a system in the real world, namely 

L = i J d'^x^''^ {x)ip{x) — J d^xd^x'il)\x)t{x,x')il){x') 

j d^xd^x'^\x)ij\x')U{\-^ - x'|)^(x')^(x), (12) 

where if) is taken from now on to represent a fermionic field. Here we shall assume the 
system is translationally invariant, so in 3- momentum space this becomes 

L = iY, ajs - Yl ^A^q - ^ (13) 
q q 

where the a, are the Fourier transforms of the field operators, depending upon time and 
momentum, and with U we have denoted the Fourier transform of the potential part of 
the lagrangian; is usually taken to be the single particle kinetic energy. Two obvious 
simplifications could be to approximate this single particle energy either as zero or as a 
step function: attractive below the Fermi sea and repulsive above. These two cases are 
intrinsically different, as we shall see in detail in the following. We can handle them together 
by setting 

M{q - kp) - AOikp - q) ; (14) 

the former case is obtained from ([T^ ) by taking the limit A ^ 0. 

These two possible choices deserve further comment before proceeding. The former has, 
as we shall show, the great advantage of being able to reproduce the FHNC diagrammatic 
scheme. The latter instead will provide the same diagrams although with different coeffi- 
cients, but nevertheless displays two advantages: namely, it provides a naturally well-defined 
ground state and furthermore the calculations, when dynamical pions will be introduced, 
will turn out to be more appropriate. This is because the latter case allows for a more 
general choice of the two body term U. 

For the sake of definiteness we write the lagrangian in the form 
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(15) 




and its corresponding euclidean version 



Le = Lit 



it) , 



(16) 



and then we identify Le with ho- We mention here that the calculations, because of subtleties 
which we shall detail in the next section, are to be performed at finite (3 and nonzero A and 
that only at the end of the calculation shall we decide in what order we want to take the 
limits. 

To summarize, the fundamental idea in the present scheme (and the same happens in 
the FHNC) is that the drastic assumption made on the single particle part of Hq greatly 
simplifies the diagrammatic expansion of @ because all frequency integrations become triv- 
ial and the four-dimensional problem reduces to a 3-dimensional one with constant energy 
denominators. 

Finally, we should mention that once we choose an ho that is sufficiently complicated 
— and the presence of a two-body term ensures this — then we are forced to make ap- 
proximations in evalutating the variational energy Eq. Of course, an approximation to Eo 
need not be an upper bound on the true ground state energy; and thus we need to have 
a sufficiently robust method for summing diagrams in an expansion in ho to allow us to 
control the approximations made and to ensure that the calculated energies converge sys- 
tematically toward the true variational energy Eo- The FHNC scheme has been successful 
in this respect, at least in the case of spin- independent correlations, and so we model our 
approach upon that scheme. 



The average value of a local, two-body potential V taken over a Jastrow correlated wave 
function can be written as 



IV. CONNECTION WITH THE FHNC EXPANSION 




(17) 
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where is the Jastrow-correlated approximation to the ground state, whose wave function 
is 

. . . , Tat) = G(t)Q{ri, . . . , Tat) . (18) 
Here G = e^/^ is the Jastrow correlation operator 

^ = 11/^' (19) 

\(f)o) is the Free Fermi Gas ground state, and (/)o(ri, . . . , tat) is the corresponding Slater 
determinant wave function. 

We wish to show that the same result comes out from the formula 

{V) = hm i / D[^\^P] (20) 
J d^xd^Vix - y)ij\x)ij\y)^{y)ijix)e~-^-!^'^'^"^''^"^ 

where V{x — y) = V(|x — y\)5{xQ — yo), provided a suitable connection between U and JJ is 
found. At first sight the connection is simple. Since V is time independent, we can replace 
(0) with 

(V) = lim / D[iIj\iP] f d^xd^yV{\x-y\) (21) 
^t(x, 0)^t(y^ o)V^(y, 0)^(x, 0)e-G%'^'^-^^\ 



which can be rewritten, using the definition of the path integral, as 

(0o|T (v{0)e'^-i'f- "-"'"^ |0o) = (0o|e-/o'"'^-'^-V^e-^-V2^-'^^|0o) (22) 
and, since U is also time independent, (0) coincides with (|1^ if we take Le = —U/j3, that 



U = -Ua and A = . (23) 

Note that we have not invoked any commutation properties between U and V , and thus the 
expressions (|l^ and (^0|) are also equivalent if we replace V by the kinetic energy T. These 
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are in fact particular instances of a more general equivalence. For example, we could replace 
V in these equations by a nonlocal two-body operator 

O = Okpk'p'alalapfak', (24) 

kpk'p' 

and the expressions would still be equivalent. This implies, by means of functional derivation 
with respect to O, that the two-body density matrices derived from (|l^ and ( pOD also 
coincide. It follows that the averages of the kinetic energy from the Jastrow ansatz and 
from the path integral approach also coincide. We limit our discussion here to a local, 
two-body operator for simplicity. 

In writing down (^) we have chosen a particular set of boundary conditions for our 
path integral. This shall be discussed presently, but it is worthwhile here to first remind 
the reader of our goal. The point is that ( |2D| ) is much more general than (p!7|); in particular. 



there is no need to restrict the potentials \^ or f/ in (|20D to be independent of time, whereas 
this is necessary if we use V and U in (just as we call ho = Le the trial hamiltonian, we 
shall generally refer to the correlations U and U as potentials). Thus our ultimate goal is to 
develop a diagrammatic scheme for calculating (pOl) . Here we are concerned with showing 
that in the particular case when we choose a V and a U that can be handled by the FHNC 
scheme arising from (0), then for the choice (^) our new diagrammatic scheme will give 
exactly that of the FHNC. 

It is natural in this respect, starting from (pO]), to introduce a perturbation expansion in 
U. First we define the generating functional in the presence of external sources 

Z[v,r^^] = [ Z)[^t^^]e^r;;jM*)+'^n*M*)+'^^W^W]'^* (25) 



and its euclidean counterpart: 

ZEhv^] = Z[r],r]^] l^^^, (26) 

so that the average of the potential V reads 

{V) = (27) 
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5n{x) 5v(v) 6vUv) 6vUx) 



6ri{x) 6ri{y) 67]'^ (y) 6ri'f{x) 

Note that whether one uses a quantum-mechanical (real time) or a statistical (imaginary 
time) approach is irrelevant in determining a static quantity (a real number) like the average 
of V. Thus we can use indifferently the generating functional or the partition function in 
the following, at least as long as we are interested in static (ground state) properties of the 
system. The study of linear responses should instead be dealt with only in the quantum- 
mechanical frame. 

The perturbation expansion is now carried out on the generating functional by means of 
the relation 

where U{x — y) = f/(|x — y\)6{xo — yo), the "unperturbed" lagrangian is Lq = L + U, with 
L given in (0), and 

Zohv^] = [ D[V.t,^]e^/-.y"°W+^^W''W+''^W^«]'^* . (29) 



The reason we have repeated here these straightforward manipulations is that they point 
out two subtleties of the procedure. 

The first subtlety is that the boundary conditions are exploited in explicitly evaluating 
the path integral in (^), that is, in defining the inverse of the lagrangian Lq. This free 
propagator in momentum space reads 

y^[uJn, q, P) = . , . H — r (30) 

cUn - A + zr] cUrt + A - zry 

with discrete frequencies because the time integration is limited: 

'^n = y • (31) 

Here we wish to choose antiperiodic boundary conditions in order that the path integral is 
equivalent to a trace over fermionic states, and so n is restricted to be an odd integer. 
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The second subtlety, strictly related to the first, concerns the degeneracy of the unper- 
turbed ground state: as long as A is kept finite an unique ground state exists and coincides 
with the Fermi sphere, and the definition ( pOf ) makes sense both for finite (3 and in the limit 
j3 —>■ oo. Then a perturbation theory can be performed by means of eq. (pHf). If we put 
instead A = from the very beginning, any independent particle state is degenerate with 
respect to Lq and fixing the ground state becomes arbitrary. Thus, taking the limit f3 oo 
for A 7^ picks out the ground state when the boundary conditions on the path integral are 
chosen as above, that is, so that (^9|) gives a trace over states. For A = 0, however, we must 
fix the ground state "by hand" — that is, as is detailed in Appendix B, we must choose a 
more involved set of boundary conditions which explicitly give the states (pQ at t = —{3/2 
and 0Q at t = /3/2 in (pH]). Moreover, with A = the momentum space propagator (|30D is 
then singular in the (3 ^ oo limit. 

To demonstrate that we are truly evaluating the average of V as given in (|T^ we need 
the expression of the propagator in configuration space, namely: 

Gl{x,y-(3) = (0o|T(^(x)V^t(^)) 10^) (32) 

_ p-iA\xo-yo\ /■ „ik-(x-y) 

J {27ry 

[e{xo - yo)e{k - kp) - e{yo - xo)9{kF - k)] . 

This corresponds to the second choice of boundary conditions discussed above, so that we 
can properly take the A — > limit. Thus the Fourier transform of ( p2D is equal to (0) only 
in the (3 oo limit (which must be taken in the end, in any case). 

With some elementary manipulations and defining the matrix element of the free density 
matrix in configuration space as 



p(A;^|x - y|) = (x|p|y) = / j^,e^'"^''-MkF - |k|) , (33) 
we can rewrite ( ^2] ) in the form 

Gl{x, y- /3) = e-^^l^»-^»l {e{x, - yo)5(x - y) - p(A;^|x - y|)} . (34) 
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In the euclidean world, one has 

6;° ^(x, y- /?) = e-^\^^-y°\ {e{x, - y,)5{^ - y) - p(fc^|x - y |)} . (35) 

To complete the proof that this perturbative expansion leads to the FHNC diagrams it is 
now sufficient to follow the derivation of the FHNC diagrams given inlll: there it is shown 
that the FHNC diagrams are obtained from a Feynman diagram expansion followed by time 
integration (we shall show this in some detail in Sec. 6). The identification (pSf) provides the 
link between the potentials, and the free Green's function is nothing but ^A£;(^5Z/j 1^) Ia=o- 
Thus we can conclude that the FHNC diagrammatic scheme exactly follows from the 
perturbation expansion (pHj) in the limit A followed by the limit /? — > oo. We stress 
that it need not be equivalent — and this is indeed the case — if one exchanges the order 
of these limits. We can infer the reason for this trouble. It is apparent in fact from (|35D 
that two cutoffs are present, the first being (3 and the second A, and dropping the first or 
the second need not necessarily give the same results. It will be proved in the next section 
that the two expansions obtained by taking the limits in different orders are in fact not 
equivalent. 

V. THE ORDER OF THE LIMITS (3 ^ oo AND A ^ 

The fundamental task in constructing the FHNC diagrams is to perform, in some way, 
the time integrations in the Feynman diagrams, in order to get a static or 3-dimensional 
theory. As a first step we prove, as a byproduct, that the two limits we discussed before are 
indeed not equivalent. We consider for simplicity in this section only vacuum to vacuum 
diagrams, but of course the procedure can be easily generalized to any kind of Green's 
function. 

To start with, let us consider a generic Feynman diagram of order n (that is, a diagram 
that contains n potential lines). It will contain, of course, 2n fermionic lines which we 
write in the form (|52D, in which the different time orderings are explicit. We shall work in 
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configuration space, and we are now interested in carrying out the time integrations. The 
general form of the diagram will be 

^ 2 2 

- j d^Xi. . .£xn j dTi. . . j rfr„ J|f/(|xfc -x;|), (36) 

2 2 

where we have divided by 13 in order to take the (3 oc limit. 

As a first step we can break the time integration region in n! parts, each one characterized 
by a definite ordering of the times. For a given ordering in eq. (|36|) , then, only one of the 
two terms composing the Green's function (^) can contribute. So each line is characterized 
by its propagation forward or backward in time, and the n! terms so obtained are, in general, 
different — although, in practice, many are equivalent. We will refer to them — improperly 
— as Goldstone diagrams. We observe that all Goldstone diagrams derived from the same 
Feynman diagram have the same spatial part, and it is only the time integrations which are, 
in principle, different. 

Having chosen a Feynman diagram and having expanded it in Goldstone diagrams, let 
us further select one of these diagrams and perform its time integrals. The diagram is 
characterized by n time variables, in a given time ordering. Let them be ri < r2 < . . . < r^. 
The ^-functions will then be absorbed into the integration limits. The time attenuation 
factors e"'^''^'^'^-'' remain, but now all absolute values can be removed. It can be easily seen 
that in any diagram, for each potential line taken at a time r, each propagator connected 
to the potential line from below (no matter if incoming or outgoing) carries a factor e~^'^ 
and each propagator connected from above carries a factor e^"^. Now let us make a trick. 
For each fermion line crossing but not stopping at a given time Tj let us insert the factors 
gArig-ATi_ 'j]2g]2 let us collect together all the factors e~^('^»+i"'^»). Their number is just the 
number of fermion lines flowing between Tj and Xj+i (including, then, also those which do 
not end at Tj or Xj+i). Thus let us call ni the number of fermion lines flowing between ri 
and T2 and, more generally, those between Xj and Xj+i. The time part of the diagram will 
then be 
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2 Tn T3 T2 



liP, A) = ^ j dTnJ dTn-1 ... J dT2 J dnll e-"«^(-«+i-0 . (37) 

^ n R n n i=l 



2 2 2 2 



Now consider the two different limits: 
1. A ^ 

In this case the integral is trivial: we get 



I{P,A) = ^. (38) 



2. p^oo 

Changing variables to Aj = — Tj+i and r„ — > /?r„, we find 



„„i 

hm A) = / dr^ f rfA„_i ... / d\^\{ e"'^^' = IT — T • (39) 

-1/2 -oo — oo ^ 

We are now ready to prove the inequivalence of the two limits. 

Consider first the diagram of fig. |^: it generates 6 Goldstone diagrams, all equal and all 
having the same graphical representation as in fig. |l]. The time integrations are immediate 
according to our preceding rules and, accounting for the multiplicity of the diagram, they 
provide a factor for the limit A ^ and for the limit (3 ^ oo. 

Next consider the diagram of fig. 0. It generates 24 Goldstone diagrams: this time, 
however, not all are equivalent. As a matter of fact, three kinds of diagrams arise, each with 
multiplicity 8, as shown in fig. ^. Their evaluation is now simple according to our rules. For 
the sum, we find (3^ in the limit A 0, of course, and 

(4A)3 ^ (4A)3 ^ (4A)2(8A) 16A3 

for the limit (3 oo. It is thus clear that no linear relationship between A and (3 can be 
found which is able to provide the same result for the two graphs. This counterexample 
proves the inequivalence between the different orderings of the two limits. 



18 



In fact, it is clear from ([36| ) that by re-scaling U ^ U/ (3, the time integrations give a 
factor I{a) = A)//3"~^, that is, a function of just the product a = (3 A. Thus we could 
send (3 —>■ CO and A ^ in such a way that a is finite, and then treat a as a new variational 
parameter. One can see from (^) that the time integrations are then in general nontrivial. 
Of course, the results P5| ) and ( P^ ) correspond to a — and a — oo, respectively; and, 
because the time integrations simplify greatly in these limits, these are the only cases we 
shall consider in the following. 



VI. DERIVATION OF THE FHNC DIAGRAMS 



To carry out the time integrations and to derive the FHNC diagrams we start from the 



separation (0). It is convenient to introduce new symbols in the diagrams. Thus let us set 

KHx; (3) = e-^^l^°le(a;o)<5(x) , (40) 

graphically denoted from now on with a double solid line, and 

^°(x;/?) = -e-^^l^°lp(fc^|x|), (41) 

denoted by a single line, so that = ]C% + g%. The procedure we shall follow is now, from a 
given Feynman diagram, to expand it into products of /C^ and g%, thus generating a certain 
number of new diagrams — which we shall call "intermediate diagrams" — with single and 
double lines, as exemplified in fig. ^. Then we shall re-sum classes of these new diagrams 
in order to simplify the subsequent time integrations. In fact, for the case A = 0, we can 
sum together the intermediate diagrams in such a way as to eliminate all time dependence 
in the fermion lines, so that the time integrations are trivial and the resulting diagrams are 
equivalent to those of the FHNC expansion. 

In order to reproduce the FHNC, of course, we consider potentials U (corresponding to 
Jastrow correlations) which are independent of time. Thus the only time dependence in the 
diagrams we are considering comes from the fermion propagators. (For the case A=0, in 
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particular, we see from (|40| ) and pi]) that the only time dependence is in the double line 
propagators /C^.) Thus, if we can demonstrate that we can sum up intermediate diagrams 
in such a way that — without ever performing any integrations — time dependence can 
be eliminated in an arbitrary line of fermion propagators, then we will have succeeded 
in demonstrating that all the intermediate diagrams can be summed to form a new set 
of diagrams which are entirely time independent. Again, for the case A = 0, it is in 
fact sufficient to consider a fermion line consisting of an arbitrary number of double-line 
propagators preceded and followed by a single line propagator. We shall begin with the 
simplest case. 

To get rid of the time dependence of the intermediate diagrams, we want to cancel the 
^-functions carried by the double lines. To this purpose consider first a double line in a 
generic diagram preceded and followed by a single line, as in fig. Then another diagram 
differing from the first only in the subdiagram of fig. ^3 surely also exists, independent 
of what is contained in the rest of the diagram. The subdiagrams plotted in fig. ^ may 
be summed immediately to provide (in Euclidean space, and with the obvious notational 
shortcuts pij = p{kp\yii — Xj|) = pji and 6ij = 6{\xi — Xj\)) 

_^^^A\t,^ts\^~A\t,-U\Q^t^ _ i^^j _ (42) 

Two points in this equation deserve comment: first, the spatial part of the diagram is 
factorized, owing to the spatial (5-function contained in the double line propagator and, 
second, the time structure is complicated and depends explicitly upon the time ordering, 
as one can easily verify. It clearly simplifies drastically in the limit A ^ 0, however, when 
it becomes just 5i2Pi3Pi4- In this limit, then, the sum of the subdiagrams in fig. 5 is 
a new, time-independent subdiagram whose spatial part is identical to that of these two 
subdiagrams. 

The next step is to consider the subdiagrams of fig. ^. It is easily realized here that, 
always in the limit A — > 0, again all the spatial parts factor out and the time part becomes 
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a sum of ^-functions, just the sum of all possible time orderings, and it gives exactly 1. 

The procedure can be iterated in this way for any number of subsequent double lines. 

Other particular cases are straightforward. For instance, if the two external lines join at 
the same point (i.e., = t^) clearly nothing changes, and the proof follows the same path 
if the external lines are replaced by a single external line closing the loop (it is sufficient 
to multiply the diagram by [g^] ^ (x^ — X3) and to integrate over ^3 and t^). Since any 
closed loop must contain at least one single line (because only the single line pieces carry 
an advanced part) we have proved that, collecting a suitable class of Feynman diagrams 
together, in any closed loop the ^-functions of /C^ cancel out. Since now any closed diagram 
contains closed loops only, our proof holds true for all closed diagrams. In this case all time 
dependence disappears, and the time integrations are trivial, giving a factor /?"■. Note that 
this proof does not depend upon the particular form of the dashed (potential) lines. They 
represent, in general, the two-body part of the trial hamiltonian ho, namely U = —U/f3 (the 
factors ^ cancel the /3's coming from the time integrations), but any one line in a given 
diagram should be replaced by V, namely, the true potential of the system. In this way the 
closed diagrams just represent the average value of V in terms of time-independent diagrams. 

Before concluding the discussion of the diagrammatic scheme in the limit A — ^ 0, we 
want to stress a last detail. It is not always straightforward to determine which diagrams 
must be collected together in order that the ^-functions sum to 1. The needed diagrams 
surely exist, but they could also coincide. That is, a given symmetry could mean that n\ of 
the necessary diagrams might not correspond to distinct Feynman diagrams. In this case 
we simply get a factor l/n\ in front of the diagram. This is in particular the case for the 
ladder diagrams of fig. |^. There, clearly, accounting for the factor 1/n! and for the locality 
of the propagators /C^, the sum reads 

W(|x - y|) + iw^dx - y|) + iw^'dx - y|) + . . . = - 1 . (43) 

This provides the still unspecified link between U = —-^U and the variational parameter 
f{rij) used in the FHNC There, in fact, no ladder diagrams are present, but those of fig. 7 
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can be summed up to give the result (^31), so that the desired hnk is 

/2 _ 1 = e" - 1 . (44) 

Having thus summed together the various classes of intermediate diagrams to obtain 
a new set of time-independent diagrams, we may shrink the double-line propagators — 
which now just represent spatial (5-functions — to points. We may furthermore replace each 
potential line with the ladder sum — 1 above. Although more than one potential line can 
meet at a single point, we must now specify that any two given points can be connected by at 
most only one potential line so as to avoid overcounting. There is one spatial propagator pij 
entering and one leaving each point (including the possibility pa), and all diagrams contain 
only closed fermion loops. This completes the derivation of the FHNC diagrammatic scheme. 

Let us now make the following observations. First, it is clear that for A 7^ our re- 
summing of intermediate diagrams does not eliminate all time dependence, and we will 
need to look for a modification of FHNC techniques in order to handle the nontrivial time 
integrations in this case. Second, we should like to stress once again that our derivation of 
the FHNC diagrams in the case A = is entirely independent of the form of the potentials. 
In particular, the unlabelled ends of the potential lines in figs. 5 and 6 could be fixed at 
different times without disrupting the conclusions that the fermion lines in the re-summed 
intermediate diagrams are time independent. Thus our approach allows for the possibility 
that the true potential V and/or the correlations U can be time dependent. Because the 
fermion lines are nevertheless time independent in this latter case, the time integrations of 
an nth order diagram become simply independent averages of the n potentials over time. 
Consequently, for the case A = 0, there will be no gain in allowing the correlations to depend 
on time; nevertheless, our formalism now allows us to find variational upper bounds on the 
ground state energies of time-dependent effective hamiltonians, and no changes in the FHNC 
machinery are needed to implement this extension. We shall see a specific example of this 
latter case in sec. 8. 
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VII. THE DIAGRAMS FOR FINITE A 



Consider now the case of the hmit (3 ^ oo followed by A — 0. This time let us make 
the rescaling U AU and consider any diagram of order U^~^ in the expansion of the 
average of V. After the frequency integrations in momentum space, it is easy to see that 
there will be iV — 1 energy denominators which are just multiples of A (this can also be 
seen from our configuration space analysis in sec. 4). Since the diagram also carries — 1 
U lines, after the rescaling all the A's cancel so that the result is independent of A. Thus 
the limit A — is unnecessary. This is, in fact, just a re-statement of our conclusion from 
sec. 5: namely, that it is only the product a = Af3 which is important here. The FHNC 
diagrammology is reproduced by taking the a ^ limit, and the present case corresponds 
to the limit a oo. 

Apart from the different rescaling, the derivation of the diagrams follows precisely the 
same steps as in the previous section but, having separated the /C^ and the g% lines, the 
frequency integration provides, as seen before, different coefficients. Once these have been 
found by explicit integration, then all IC^ lines disappear and, up to the different coefficients, 
we are left with exactly the same diagrams as before. This time, however, the re-summation 
of the class of diagrams of fig. follows a different path. Let us work in momentum space 
and let us define U{q) as the Fourier transformed potential and T{p, k; q) as the sum of the 
ladder diagrams, p and k being the momenta entering the diagrams. Clearly r(p, k] q) fulfills 
the Bethe-Salpeter equation 



Going to configuration space, this equation is immediately solved. Calling xi, X2 and yi —1/2 
the transforms of p, k and q, we get 



T{p,k;q) = U{q) + i 



I 



dH f{p,k;t)U{q-t) 



(27r)4 (po + to - A + iri) (A;o - to - A + irj) 




(45) 




(46) 



—00 



—00 
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^ f/(yi-y2) 



h + po - 2A + if] 

which still depends upon the initial times. 

However, we recall that no conditions have been imposed upon the possible time (or 
energy) dependence ofU. In particular, we can rewrite the trial hamiltonian (now, strictly, 
no longer a hamiltonian) as 

ho= J2 ^4«fc - J2 A^Ofc + A^A(/c,p;q)aj4afc_qap+g (47) 

k>kp k<kp kpq 

where we explicitly allow an energy dependence in the "effective potential" A upon the initial 
energies and po but not upon the transferred energy go- We can now write down the Bethe- 
Salpeter equation exactly as before, getting a similar result. Going again to configuration 
space and calhng A the Fourier transform of A with respect to the spatial coordinates only, 

we get 

r(x„ x„ y„ y2)= J^l ^e-(^'?-?)-^^^o(.r..o) (43) 
—00 —00 

A(Po,fco,y,-y,) _ _ 

^ A(po, ko, yi - ya) 
^0 + Po - 2 A + ir] 

Now we can use the extra energy dependence of A to remove that of T. We can define, in 
fact, 

w , X jpo + h- 2A) 7(y) ..Qx 

A Po, ko, y = — — — — — (49) 

Po + ko-2/\ + 7(y) 

so that inserting (^) into (^T]) we get the straightforward result 

r(xi, X2, yi, 2/2) = 7(yi - y2)6{xl - yl)5{xl - y'MVi " v') ■ (50) 

Thus, having chosen a particular frequency dependence of the effective potential A, we can 
substitute everywhere the potential with the sum of the ladder series F, which is now fully 
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frequency independent. Of course at this stage the only relevant point is the existence of 
the starting quantity A, not its cumbersome form. Naturally, the contact with the FHNC 
diagrammatic scheme is now provided by — 1 = 7. 

However, we must remind the reader that in the (3 00 limit in which we are working 
now, the coefficients of our diagrams resulting from frequency (or time) integration are not 
those of the FHNC because the fermion lines are not time independent. We still have to 
make the connection with FHNC techniques for summing diagrams in this case, for the 
different coefficients prevent us from summing diagrams through the usual FHNC integral 
equations. While we shall not pursue this topic further in the present work, we will mention 
that this A 7^ case might be handled by extending the FHNC equations to include time 
integrals, thereby generating the different coefficients as one builds chains and hyperchains. 
It is convenient to work in configuration space, but to transform from time to frequency 
variables. Then one arrives at a diagrammatic scheme which is identical to that of the 
FHNC except that exchange and dynamical correlation lines carry frequency variables, and 
the points — while still representing spatial 5-functions — now have a nontrivial frequency 
dependence. Because of the simple form of the single particle term in Lq, however, one can 
still classify diagrammatic types in a way similar to that done in the FHNC, forming chains 
by performing convolution integrals over frequency as well as spatial variables, and forming 
hyperchains by summing over ladder diagrams as demonstrated here for the simplest type 
of two-point correlation. An important extention here would be to allow for correlations A 
in (|47|) which also depend on the frequency go- 



Finally we come to a quantum field theoretical model in order to see how it can be 
embedded in the previous formalism. Consider for the sake of simplicity a system of non- 
relativistic nucleons interacting with a pionic field through the typical interaction lagrangian 



VIII. THE PION EXCHANGE POTENTIAL 




(51) 
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If Lat and L,r are the free lagrangians for nucleons and pions, respectively, the partition 
function of the system reads 

D[V'V0]e"^-^/2''"^'''^^^"^'"^*^'^ (52) 

where it must be stressed — and this is the main drawback of the present model — that 
the term in the lagrangian has been neglected, whereas we know that it is necessarily 
present in a renormalizable theory. This simplification ensures, however, that the remain- 
ing lagrangian is quadratic in the pion field, which can consequently be integrated out. 
Calculations are straightforward^ and provide (apart from an inessential constant factor) 

Dltfj^iP] (53) 
where the fermion current is 

iUx) = -t^\x){(TT'')^ix) (54) 

and Do{x — y) denotes the free pion propagator, which in momentum space and in the 
euclidean world reads 

ujl + q^ + ml 

with LVn given of course by (pi]), this time with n even integer only. We can then follow all 
our machinery to derive the average value of the interaction term, which looks exactly as 
before if we disregard its frequency (or time) dependence. 

Now from (pTj) we observe that the diagrams we need to evaluate contain one pion 
propagator (hence a line with the extrema put at different times) closed by a two-particle 
Green's function whose time evolution is governed by the trial hamiltonian. Let us now 
take the limit A ^ 0, j3 ^ oo {a 0) as in sec. 6. We recall that the rules for obtaining 
the FHNC diagrams derived there are not influenced by a possible time dependence in the 
interaction lines. More explicitly, it is irrelevant if the potential lines in figs. |^ and |^ are 
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instantaneous or not. In both cases, for a given Feynman diagram, others exist (including, 
perhaps, equivalent diagrams) such that when these are added together the 6'-functions 
sum up to 1, and all time dependence drops out. Consequently the two-particle Green's 
function evaluated in the frame of FHNC turns out to be time- independent. Let this be, 
say, G//(xi, X2, X3, X4). The integral we are interested in is then 
1 r rl^/'i /•/3/2 

-Tr / Sxd^y / dxo / dy^ (56) 
De{x - y) {-ia- ■ V^t) {-i(T ■ Vyr) Gii{x, y, x, y) . 

Now we carry out the integrations over the imaginary times. Using eq. (^), it is a simple 
matter to show that 

1 /•/3/2 |-/3/2 

- / dxo dyoDsix - y) (57) 

" y _-jq (x-y) ^ 



/•/3/2 


/•/3/2 r 


/ dxQ 


/ dyQ / 


1-13/2 


'-/3/2 J 



re 



E 



" ^ ^-^q-(x-y) ^ _ ^ 



re 



(27r)'^ + ml 47r|x — y| 

The conclusion is that what we need for the calculation is to take the effective potential at 
zero energy and then to perform the standard FHNC calculations. In the present model, 
the effective potential is just the OPEP, as one might have expected. 

To go beyond this variational approximation, now, we have several possibilities. We 
can try to improve the variational result by allowing the trial hamiltonian to have a spin- 
dependent interaction and/or a time dependence by taking A 7^ 0. Otherwise, we can keep 
the simpler trial hamiltonian used here, which allows the use of the full FHNC machinery, 
and we can calculate corrections to the variational result through perturbation thery as 
described in sec. 2. Implementation of these improvements is currently underway. 



IX. CONCLUSIONS AND OUTLOOK 



Let us now briefly summarize the results obtained in this paper and discuss the possible 
developments they naturally suggest. 
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The central idea in this paper is that of formulating a variational approach for dealing 
with strongly interacting systems which goes beyond the limitations of the Fermi Hyper- 
netted Chain and Correlated Basis Function methods. This led us to use a path integral 
approach, where the variational method consists in choosing a trial hamiltonian instead of 
choosing an ansatz for a nonrelativistic wave function. We have obtained several important 
results. First, one can define the path integral and choose the trial hamiltonian in such a 
way that for particular systems for which a nonrelativistic wave function makes sense, our 
approach produces exactly the same results as the successful FHNC formulation. Second, 
in this path integral formalism one is not linked to a true hamiltonian, but an effective 
hamiltonian can be used as well. Thus ultimately we can generahze the FHNC variational 
techniques for use in studying the enormous variety of strongly correlated systems for which 
the usual FHNC cannot be directly applied because a formulation in terms of nonrelativistic 
wave functions breaks down. Thirdly, we have shown how corrections to a static variational 
approximation can be added by choosing a more general trial hamiltonian so to improve the 
variational energy and by including perturbative corrections in the difference between the 
true and trial hamiltonians within a compact path integral formulation. 

Having first briefly described the variational principle in the Feynman path integral for- 
mulation, we proceeded to define our approach by using the FHNC and CBF approaches 
as motivation. Since these latter approaches treat correlations between particles as more 
important than independent particle kinetic energies, in our approach the fundamental sim- 
plification with respect to the true hamiltonian h we chose for our trial hamiltonian ho was 
to modify its one-body piece. For much of the paper we limited ourselves to the case of a 
time-independent ho, and throughout we only considered spin-independent interactions in 
ho- We made a specific simple choice of the one-body term in ho which depends on a param- 
eter a — /3A, and we have shown that in the particular case a — > our expansion scheme 
in terms of Feynman diagrams exactly reproduces the FHNC results for a variational wave 
function of the Jastrow type. We also showed that in the limit a ^ oo, the time integrals 
are also simple, if not as trivial as in the case a — 0. We showed that this other limit 
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corresponds to a new approximation scheme which has the same diagrammatic structure as 
the FHNC except with different coefficients in front of each diagram. 

Finally, to show how our method works in practice, we have shown that the variational 
scheme can be applied to a quantum-field-theoretical case, a system of interacting nucleons 
and pions. We have proved a minimal but highly nontrivial result — namely, that in a 
system containing dynamical pions a variational procedure built up upon the static one 
pion exchange potential still provides an upper bound to the energy. 

The present paper was mostly devoted to proving the correctness of the method and 
its equivalence, in a very particular case, with the available FHNC calculations. However, 
it opens several interesting perspectives which we want to discuss now. First, we limited 
ourselves to trial hamiltonians which correspond to the Jastrow ansatz for the wave function 
in the usual FHNC approach. State-dependent interactions in Hq can be handled on the same 
footing as those described here. Very likely, as some preliminary results seem to indicate, 
the method presented here could provide a very useful tool to overcome the limitations of 
the Single Operator Chain approximation. 

Next, we have to exploit the new possibilities offered by the present approach. First, 
energy-dependent correlations can be accounted for, thus improving upon the usual cal- 
culations carried out with phenomenological nuclear potentials. Second, an extension to 
relativistic systems is now allowed, potentially cumbersome but in principle feasible. In 
fact, again we are faced, in a relativistic system, with a case where a trial wave function 
looses its sense but our (eventually time-dependent) trial hamiltonian does not. Third, in 
connection with the previous point, the possibility of carrying out a variational calculation 
with the Bonn potentialEl is now offered. Note that using energy-dependent interactions in 
the trial hamiltonian (energy-dependent correlations in the FHNC language) allows for an 
improvement upon the results of sec. 8, going beyond the static limit of the potential. 

Finally, another promising application of our approach is to quantum field theory, as 
indeed we briefly demonstrated in sec. 8. Variational methods have been little used in 
field theory up to now (see, for example ,1^1-0), and to our knowledge nothing of the sort 
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introduced here to deal with strongly correlated systems has been explored in the field 
theoretical problems. Our approach thus opens up new possibilities in this direction. The 
obvious, if ambitious, application would be to study low energy bound states of quantum 
chromodynamics, where the usual perturbative methods are hopeless. There are two paths 
one can follow in this respect. One is to study euclidean path integrals as in sec. 8, obtaining 
an effective lagrangian through functional integration. A second possibility, which might 
overcome some of the problems mentioned by FeynmanS, is first to develop an effective 
hamiltonian using Wilsonian renormalization group methods — a promising approach in 
this direction is proposed and detailed in§ — and then to analyze this effective, frequency- 
dependent hamiltonian with our path integral variational methods. Indeed, it was the 
availability of powerful hamiltonian methods such as the variational principle which helped 
lead Wilson to develop his renormalization group approach^. In either of the approaches 
to field theoretical systems mentioned here, of course, one has difficult renormalization 
problems to overcome. 

It is our hope, then, that the variational approach described here will eventually provide 
an important tool for studying many different types of strongly correlated systems for which 
the usual perturbative techniques are insufficient, and that only a small number of the 
potential applications of our approach have been mentioned here. 
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APPENDIX A: DERIVATION OF EQ. (8) 

We show here that the first and third terms on the r.h.s. of eq. (0) cancel even for a 
time-dependent potential h^r). 

Let us define a suitable "partition function" 
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SO that 



(Al) 



By comparison with (|^) 



ID 



13/2 \ 
ho{T)dT) ) 

/ho 



(A2) 



A=l 



(A3) 



Now we represent ho{T) in ( [A2[ ) as the Fourier series J2uj„ ^^"^"^^hoiun) with ujn = T^nj (5, so 
that as /5 — s> oo 

f/3/2 



/3/2 



/io(r)rfr — ^ /?/io(0). 



(A4) 



In this hmit, then, the path integral is equivalent to a trace over the eigenstates of /io(0), 
whose eigenvalues we write as e„. Thus 



1 a_ 



log^o' 



/D[V^t^]/ip(0)e-^/3ho(o) 
/L)[7/;tV;]e-^/3/ro(o) 

71 "^n*-^ 



(A5) 



As /3 — ^ oo, of course, the lowest energy state dominates this average, giving the result 
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which does not depend upon A. This implies 



since Z is normalized such that log 
with 



A=0 



0. Comparing 



A=l 



and 



1 / /-/^/a ^ 
Eo< Mm -( hdr 

P^oo (3 \J-p/2 I 



(A6) 
we end up 

(A7) 



ho 
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APPENDIX B: BOUNDARY CONDITIONS ON THE PATH INTEGRAL 



Here we show explicitly the two choices of boundary conditions on the ferm-ionic eu- 
clidean path integral discussed in Sec. 4. Further details may be found, for example, in 
Ref.0. 

The standard notation 

D[^^^]e-^-^/^'°'^'''^^'" (Bl) 
used in the text is a misleading shorthand for 

/ f/3/2 
lim n d^\z{i5nnN)\e~^-^i^'"'^"'"*^''\ (B2) 

z{-f3/2)=z, "=-^ 

where the integration measure is 

dij[z] = Yldzldzqe-^'". (B3) 

Here ho is assumed time-independent for simplicity, and the classical Grassman variables z 
and z* correspond to the nonorthogonal coherent states 

1^) = e^9^'"^|0) and (^1 = (0|e^«^'"% (B4) 

with |0) the Fock vacuum. These coherent states are nonorthogonal superpositions of Fock 
states with differing numbers of particles, and their overlap is 

= e^«''"' =e"*"'. (B5) 



The variables Zi and are not integrated over in ( |B2D , which is equivalent to the matrix 
element 

{zf\e-'^'''\zi). (B6) 
The usual choice of boundary conditions is to set z*jr = —z* and to integrate over Zi 

diJ,{zi)Z{-z*,Zi), (B7) 
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so that the generating function Z becomes a trace over all Fock states — that is, the partition 
function. Then the /3 — > cxd limit will pick out the ground state, provided has a unique 
ground state. If this is not the for the FHNC choice A = 0, we must work a little 

harder to ensure that ( |B2[ ) is equivalent to the desired matrix element. We can use the 
closure relation 

^ d[i\z\\z)(^z\ = \ (B8) 

to represent any state |0) as 

10)= ^ dix\z\\^){z\<\>) . (B9) 
Thus in the A = case we must take our shorthand expression (|B1| ) to represent 

Z = y" d\i\zj\dii\z^{^Q\z^)Z{z),zi){zi\^Q) (BIO) 

in order to ensure that we are evaluating the matrix element (0o|e~'^'^'' |0o)- Here |0o) is the 
free Fermi gas ground state as in Sec. 4, so that 

(0oi^) = n (Bii) 

q<kp 
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FIGURES 
FIG. 1. Third order Feynman diagram 

FIG. 2. Fourth order Feynman diagram 

FIG. 3. Time-ordered "Goldstone" diagrams corresponding to fig. 2 

FIG. 4. Example of separation of propagators into two terms 

FIG. 5. Two subdiagrams which sum to give a time-independent diagrammatic structure 

FIG. 6. A further example of summation to form time-independent diagrams 

FIG. 7. Summation of ladder diagrams 
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